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MHD WAVES IN SUNSPOTS 


Robert Sych^ 

1. INTRODUCTION 

Studying magnetohydrodynamic (MHD) waves 
and oscillations in the solar atmosphere is one of 
the fastest developing fields in solar physics, and 
lies in the mainstream of using solar instrumen¬ 
tation data. Wave processes draw our attention 
due to their possible role in problems of coronal 
physics, among which: coronal heating, flare en¬ 
ergy release and coronal mass ejections, solar wind 
acceleration. Because the observed properties of 
waves include the information about propagation 
medium, they are natural probes providing us with 
its diagnostics means. 

Analyzing quasi-periodic oscillations (QPOs) of 
the emission above sunspot, first discovered by 
[Beckers and Tallant, 1969] and referred as ”um- 
bral flashes” takes a special place in studying 
wave processes. Although 3-min sunspot oscilla¬ 
tions have been observed for decades with sufficient 
progress in understanding this phenomenon [Bog¬ 
dan and Judge, 2006; Thomas and Weiss, 2008], 
we are still far from realizing many problems con¬ 
cerning physics of sunspot wave processes. For ex¬ 
ample, the relation between 3-min umbral waves 
and 5-min penumbral traveling waves has not been 
clear yet. Different studies result in mutually ex¬ 
clusive conclusions: from asserting that there exist 
prerequisites for penumbral wave emergence from 
sunspot umbral oscillations [Christopoulou et al., 
2001] to addressing them as independent phe¬ 
nomena [Christopoulou et al, 2000; Kohanov and 
Makarchik, 2004]. Recent studies [Tziotziou et al., 
2006, 2007; Thomas and Weiss, 2008] showed that 
umbral flares and penumbral traveling waves are a 
manifestation of a common wave process with wave 
fronts propagating through the sunspot umbra and 
penumbra. 

It is generally accepted that magnetoacoustic 
waves move at a local sound speed along magnetic 
field lines from the photospheric layers through the 
chromosphere into the corona. However, it is still 
unclear how high they propagate into the corona. 
What determines their properties in the corona? 
How does the sunspot wave cut-off frequency vary? 
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How are the oscillation parameters spatially dis¬ 
tributed depending on the oscillation period? Are 
they directly related to propagating longitudinal 
waves usually observed in ultraviolet in coronal 
magnetic structures like EUV fans [De Moortel and 
Nakariakov, 2012]? What is the nature of the re¬ 
cently discovered link between 3-min sunspot oscil¬ 
lations and the similar-periodicity flare oscillations 
[Sych et al., 2009]? An insight into these problems 
will give us a chance to diagnose plasma at differ¬ 
ent levels of the solar atmosphere, because waves 
carry the information on the nature of their emer¬ 
gence and on the properties of the medium where 
they propagate. 

The review on dynamic of sunspots is ar¬ 
ranged as follows: Section 2 addresses the spatial- 
frequency morphology of sources of sunspot os¬ 
cillations and waves, including their localization, 
size, oscillation periods, height localization with 
the mechanism of cut-off frequency that forms the 
observed emission variability. Section 3 presents 
a review dynamic of sunspot wave processes, pro¬ 
vides the information about the structure of wave 
fronts and their time variations, investigates the 
oscillation frequency transformation depending on 
the wave energy. Section 4 addresses the initializ¬ 
ing solar flares caused by trigger agents like magne¬ 
toacoustic waves, accelerated particle beams, and 
shocks. Special attention is paid to the relation be¬ 
tween the flare reconnection periodic initialization 
and the dynamics of sunspot slow magnetoacoustic 
waves. Section 5 presents a short review of theo¬ 
retical models of sunspot oscillations. 

2. SPATIAL STRUCTURE OF OSCILLATION 
SOURCES ABOVE SUNSPOTS 

Magnetoacoustic gravity (MAG) waves above 
sunspot atmosphere have a strong dispersion when 
propagating from the subphotospheric layers into 
the corona. Their properties dramatically depend 
on the cut-off frequency that, in turn, depends on 
the propagation medium plasma parameters. The 
waves with the frequency under the cut-off have a 
strong damping, and can not reach the upper lay¬ 
ers of the solar atmosphere. Due to that, there 
occurs a filtration of the propagating broadband 
pulses by frequencies and height [Suematsu et al, 
1982; Botha et al., 2011; Yuan et al, 2013]. The 
cut-off frequency in the stratified isothermal at¬ 
mosphere with magnetic field was obtained ana- 
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lytically in [Bel and Leroy, 1977; Zhugzhda and 
Dzhalilov, 1984], The frequency value in general 
case depends on the local plasma /3 and on the 
magnetic field inclination angle ip. In the regions 
with a high plasma /3 (photosphere and chromo¬ 
sphere of the quiet Sun), this frequency comes to 
the ordinary acoustic cut-off frequency, ~ 5.2 mHz. 
For low plasma /3 (sunspots, corona), the cut-off 
frequency is modified by the value of the inclina¬ 
tion angles of power magnetic lines along which 
waves propagate. This explain the existence of 
low frequency compression waves in the corona 
[De Moortel and Nakariakov, 2012] that are as¬ 
sumed [Fontenla et al., 1993; Jefferies et ai, 2006] 
to transport the bulk of the energy as compared 
with high-frequency waves. The experimentally- 
obtained cut-off frequency can be used as a seis- 
mological instrument to study the local magnetic 
topology in sunspots. 

In review [Borrero and Ichimoto, 2011] the 
authors give an overview of the current state- 
of-knowledge of the magnetic field structure in 
sunspots from an observational point of view. 
Briefly described of tools that are most commonly 
employed to infer the magnetic field in the sunspot 
atmosphere. Taken into account the differences of 
global and local magnetic details of sunspots, fo¬ 
cusing on the implications of the current observa¬ 
tions of energy transport mechanisms for the differ¬ 
ent sunspots models. Presented different methods 
of extrapolations of the magnetic field towards the 
corona. 

Investigations in [Bel and Leroy, 1977] showed 
that long-period waves are capable of permeat¬ 
ing through the photosphere and the chromosphere 
and of propagating upward, into the coronal re¬ 
gion. The magnetic field concentration in the 
form of pores or sunspots with a complex magnetic 
structure are natural waveguides for these waves. 
The modification in the magnetic field line incli¬ 
nation is shown to be possibly an explanation for 
the origin of 5-min waves in chromospheric spicules 
[De Pontieu et al., 2004] and coronal loops in ac¬ 
tive regions [De Moortel and Nakariakov, 2012; De 
Pontieu et al., 2005]. One also observes the per¬ 
meation of the solar global p-mode into the chro¬ 
mosphere through the inclined magnetic fields at 
peripheral portions of plages [de Wijn et al., 2009]. 
Subsequent observations of long-period oscillations 
in the corona were also explained by the existence 
of inclined magnetic channels along which waves 
propagate [Marsh et al., 2009; Yuan et al., 2011]. 

Direct observations of the magnetoacoustic cut¬ 
off in sunspots by using the spectropolarimetry re¬ 
sults were first shown in [Bloomfield et ai, 2007]. 
The authors revealed that the cut-off frequency de¬ 
pends on the change in the magnetic field in ap¬ 


proaching/approximation of a weak plasma beta 
[Bel and Leroy, 1977]. However, there was no fur¬ 
ther investigation to determine how this frequency 
is related to the inclination of field lines. It became 
clear that the observed differences in the phase and 
correlation in the time series of the intensity vari¬ 
ations obtained at different heights in the solar at¬ 
mosphere look more realistic, taking into account 
the spatial distribution of the magnetic field and 
its curvature [Bloomfield et al., 2007]. 

[Kobanov et al, 2008, 2011] revealed a peculiar¬ 
ity in sunspot oscillations as a height inversion of 
power variation. At the chromosphere level, the 
3-min oscillation power maximum spatial localiza¬ 
tion is usually over the umbra. At the photo¬ 
sphere level, these oscillations have a sharp damp¬ 
ing. The revealed regularity agrees well with the 
results of the previous investigations. In partic¬ 
ular, in [Balthasar et al, 1987], no 3-min peri¬ 
odicity components were observed. In [Lites and 
Thomas, 1985], the level of the recorded 3-min os¬ 
cillations was very low, at the noise level. By us¬ 
ing helioseismology methods [Braun and Duvall, 
1990; Nicholas et al., 2004], suppressing the level 
of photospheric oscillations in active regions was 
recorded. Analyzing the brightness fluctuations 
observed in the Hinode/SOT G-range 4305A in¬ 
dicates at a decrease in the broadband bright¬ 
ness variations in the central part of sunspots [Na- 
gashima et al., 2007]. 

To explain the effect of changes in observed 3- 
min oscillations power with height two mechanisms 
was proposed [Kobanov et al., 2011]. First of them 
connected with suppress the observed line Fe I 
6569 A at the photospheric level and decreases 
the amplitude of acoustic oscillations in deeper 
regions [Lites, 1992; Sigwarth and Mattig, 1997] 
of sunspots. According to [Balthasar and Woehl, 
1983], the depression of the umbral photosphere 
with respect to the penumbra (the Wilson effect) 
is about 700 km. This difference in height can lead 
to the decrease in the 3-min oscillation amplitude 
by a factor of 2 [Lites et al, 1998] and consis¬ 
tent with the 3-min the line-of-sight (LOS) veloc¬ 
ity power distribution obtained at the photospheric 
level [Kobanov et al., 2011]. The second mecha¬ 
nism connected with the topology of the sunspot 
magnetic field. In the umbra, the oscillations as 
a slow magnetoacoustic waves guided by the ver¬ 
tical magnetic field lines upwards and reach the 
chromosphere. At the umbra/penumbra border 
the magnetic field lines are inclined from the ver¬ 
tical direction and prevent the free propagation of 
the short period waves [Bogdan and Judge, 2006]. 
In [Finsterle et ai, 2004] pointed out that the ~ 
2.4 min oscillations are reflected by the active re¬ 
gion canopy. [Kobanov et al., 2011] proposed that 
partial reflection of the waves from the magnetic 
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dome can lead to the observed accumulation of the 
wave energy at the photospheric level outside um¬ 
bra under horizontal magnetic field. This may ex¬ 
plain the relative increase of the 3-min oscillation 
power outside the umbra in the photosphere, and 
the corresponding decrease of the chromospheric 
oscillations. 

In [McIntosh and Jefferies, 2006], the authors 
studied the sunspot narrowband signal propaga¬ 
tion time, and found good agreement with the 
study in [Bel and Leroy, 1977], both in the quiet 
Sun regions (/3 S> 1) and in sunspots (/3 <C 1). 
In [Tziotziou et al., 2006], the empirical formula 
Vpeak(7’) - 1.25I4iax(7’) [Bogdan and Judge, 2006] 
in the sunspot chromosphere was used to obtain 
the value of the magnetic field inclination. In this 
formula, Tpeak denotes the peak frequency in the 
power spectrum, and I4iax for the local cut-off fre¬ 
quency modihed by the inclined magnetic field. 
In a recent study [Reznikova et al., 2012], spec¬ 
tral and phase variability in the emission intensity 
with height in ultraviolet was investigated. The 
existence of delays between the SDO/AIA tem¬ 
perature channels was shown, which indicates the 
presence of upward propagating brightness distur¬ 
bances. The cut-off frequency height variations 
were found to agree well with the conclusions pre¬ 
sented in [Bel and Leroy, 1977]. Subsequent stud¬ 
ies [Reznikova and Shibasaki, 2012] to compare the 
cut-off frequency of the observed magnetoacoustic 
waves at 304 A (SDO/AIA) with the computed 
values of the extrapolated field showed good agree¬ 
ment with each other. 

In [Kobanov et al., 2013a], direct measurements 
of phase delays between the wave speed variations 
above sunspots and faculae were performed. Dif¬ 
ferent magnetic structures show different delay val¬ 
ues related to the wave propagation speed in the 
selected region. The obtained speeds were revealed 
to dramatically surpass the generally accepted val¬ 
ues of sonic speed in the photosphere. Thereupon, 
there arises the issue of these oscillations origin and 
of the wave mode recorded at observations. 

[Abramov-Maximov et al, 2011] compared the 
velocity oscillation modes at the photosphere- 
chromosphere level and the radio emission oscilla¬ 
tions (NoRH, 17 GHz) in the transition zone. The 
spectral composition of the wavelet spectra for op¬ 
tical and radio observations showed an about 50- 
sec time delay of signals. This indicates at the 
MHD wave propagation from below upward in a 
sunspot magnetic tube. 5-min oscillations show 
the same spectra at various heights only near the 
umbral region. 

[Kobanov et al., 2013b] analyzed sunspot os¬ 
cillations using the SDO/AIA and ground-based 
optical observation data. Low-frequency oscilla¬ 


tions are shown to form sources at the sunspot 
penumbral boundary. Their shape agrees with the 
penumbral boundaries. At the coronal level, this 
component forms the magnetic fans’ outer bound¬ 
aries. 5-min oscillations are localized at the um- 
bral/penumbral boundary. The fundamental 3- 
min component is inside the umbra. A higher- 
frequency component (~ 2.1 min) is concentrated 
in individual small areas of an umbra. The authors 
calculated the wave propagation velocities: 28 ± 15 
km/s, 26 ± 15 km/s, and 55 ± 10 km/s for the Si 
I 10827A - He I 10830A, 1700A - He H 304A, and 
He H 304A - Fe IX 17lA line pairs. Measuring 
the high rate of error indicates a lack of co-phase 
oscillations at different height levels due to their 
frequency dynamics. 

Observations of the photospheric magnetic field 
oscillations [Norton et al., 1999; Balthasar et al., 
1998; Bellot Rubio et al., 2000] showed a fine field 
fragmentation in the regions with the maximal 
power of the oscillations coinciding with sunspots. 
For example, optical observations show that the 
magnetic field 5-min oscillations [Zhugzhda et al., 
2000] are localized both in the magnetic field 
isolated tubes (pores), and at the sunspot um- 
bral/penumbral boundary. In sunspots, one also 
observes short-lived blinker-type flare events with 
a similar periodicity, 400 through 1600 sec [Har¬ 
rison et al., 1999; Brkovic et al., 2000; Lin et al., 
2006]. 

[Kobanov and Chelpanov, 2014] studied effi¬ 
ciency of the oscillation mode permeation from the 
subphotospheric levels into the corona in the fac- 
ular regions near sunspots. They used both the 
ground-based data in the Si I 10827A and He I 
10830A lines and the SDO/AIA data (Fe I 1700A, 
He H 304A, and Fe IX 17lA). The analysis of the 
spatial distribution in oscillation power at different 
heights showed that, at the coronal level (17lA), 
the emission oscillation maximum within the 11-17 
min periods range is observed. This implies that 
the low frequency oscillations dominate in coronal 
loops over plages. 5-min oscillations that prevail 
everywhere in measurements of velocities in the 
facula atmosphere lower layers are notable only in 
the detected compact parts of coronal loops. Fur¬ 
ther, [Kobanov et al., 2006] revealed the spatial 
distribution of the sunspot velocity variations in 
a chevron form. These structures reflect the exis¬ 
tence of propagating waves with a phase velocity 
of about 50-80 km/s. The comparison between the 
characteristics of the 3-min umbral waves and the 
penumbral traveling waves at the chromosphere 
level revealed no association between them. A 
possibility to explain the oscillation frequency and 
the propagation velocity depending on the umbral 
center within the interaction and transformation 
of different oscillation modes is discussed. For 
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the traveling 5-min penumbral wave, the velocity 
is shown to insignificantly vary with the distance 
from the umbral border. 

The spatial distribution of the sunspot oscil¬ 
lation frequency and power has a different de¬ 
gree of insight [Bogdan and Judge, 2006], and 
was vigorously studied by using the data obtained 
with both ground-based, and space instruments 
(SOHO, TRACE, HINODE). These studies, the 
SOHO/MDI observations in particular, showed 
that, at the photosphere level, there occurs an in¬ 
crease in the Doppler velocity and a line broad¬ 
ening, but not the continuum intensity [Jain and 
Haber, 2002]. An increase in the spectral power 
was discovered in active regions with a sufficiently 
strong magnetic field. In turn, as the field further 
strengthens, a significant decrease in oscillations is 
observed, for example, in the sunspot center. The 
observations in Ha showed that steady oscillations 
(with the period of about 2.6 min) dominate in 
the umbral region and in the inner part of sunspot 
penumbra, with the presence of traveling waves in 
the penumbra within the 5.6 min range. However, 
these waves were not recorded at the photospheric 
heights [Christopoulou et at, 2000, 2001]. The 
sunspot umbra 3-min oscillations and the penum¬ 
bral 5-min traveling waves were assumed to show 
an example of reflecting the dynamics of the same 
waves generated in the photosphere and propagat¬ 
ing along different magnetic field lines with a differ¬ 
ent inclination level [Bloomfield et at, 2007]. On 
the other hand, in [Kolobov and Kobanov, 2009], 
no association was revealed between the umbral 
3-min oscillations and the penumbral 5-min trav¬ 
eling waves at the chromosphere level. Parameters 
for the 5-min oscillation mode were obtained. The 
length of propagating waves is ~ 12-30 arcsec, or 
9-20 Mm, the phase velocity varies from 28 to 65 
km/s. 

[Christopoulou et at, 2003] showed that oscilla¬ 
tions with different frequencies within the 3-min 
period range are localized in different parts of an 
umbra at the chromosphere level. One observes 
numerous oscillation peaks in the 2-3 min spectral 
range. There are three detected harmonics in the 
spectrum: the largest is near 3 min, the second 
is near 2.6 min, and the third (about 2.2 min) is 
present for a restricted time. Within existing theo¬ 
ries, the ~3 min oscillations are considered as an ef¬ 
fect of the photospheric resonator, two other modes 
are considered within the chromospheric resonator. 

In solar atmosphere take place the formation of 
acoustic wave resonators both in active and quiet 
Sun regions [Sobotka et at, 1994; Rosenthal et at, 
2002; Srivastava et at, 2008]. Eor example the 
analysis of observational data, obtained by Hin- 
ode/EIS [Srivastava et at, 2008], show ~ 7-13 min 


intensity oscillations (He H 256.32 A, EeXI 188.23 
A, and EeXH 195.12 A) in brightest core of the 
magnetic network. The authors suggest that the 
field-free cavity regions, under the bipolar small- 
scale magnetic canopy may serve as resonators for 
the waves. The field-free cavities probably have 
granular dimensions. The MHD equations in the 
cavity and canopy regions was solved separately, 
and merge the solutions at the cavity-canopy in¬ 
terface. The dispersion relation for oscillations 
shows that the first harmonic is leaky, and that 
the oscillations may propagate upwards. Numeri¬ 
cal modeling of the wave propagation through mag¬ 
netic structures in a stratified atmosphere [Rosen¬ 
thal et at, 2002] shown, that in regions with sig¬ 
nificantly inclined field to the solar normal, refrac¬ 
tion by the rapidly increasing phase speed of the 
fast modes results in total internal reflection of the 
waves. The height of this reflecting layer will de¬ 
pend on the field geometry and wave frequency, but 
in the Sun it will certainly be a highly complicated 
surface. This varying reflection height may be re¬ 
sponsible for the intermittancy seen in the oscilla¬ 
tions. [Kuridze et at, 2009] studied the influence 
of this small-scale magnetic canopy, on the wave 
dynamics. They found that the field-free cavity re¬ 
gions, under the canopy, may trap high-frequency 
acoustic oscillations. The lower-frequency oscilla¬ 
tions may be channeled upwards in the form of 
magneto-acoustic waves. 

[Tziotziou et at, 2006] found that the 3-min os¬ 
cillation frequency in the chromosphere dramati¬ 
cally changes at the umbral/penumbral boundary 
when moving from the sunspot center. Revealed 
were small bright mottles in the sunspot center 
where the oscillations with a period less 4.1 min 
tended to decrease [Nagashima et at, 2007]. A sim¬ 
ilar effect was found on the Doppler velocity maps 
in the Ca H 8542A line [Tziotziou et at, 2007]. The 
nature of these ’’dark spots” in the spatial distri¬ 
bution of 3-min oscillations at the chromosphere 
level still is not clear. 

[Jess et at, 2013] studied the role of the mag¬ 
netic field inclination angles in forming physical 
characteristics of traveling waves in a penumbra at 
the chromosphere level. An increase in the period 
oscillation depending from the umbra center was 
found. Based on applying the cut-off frequency 
mechanism, a possibility to obtain the magnetic 
field structure above sunspot by using the data 
on dominating wave periods’ spatial distribution 
was shown. It is concluded that there is a re¬ 
lationship between the geometry of the magnetic 
field connecting the photosphere with the chromo¬ 
sphere and the behaviour of traveling waves in a 
penumbra in the upper chromosphere. This con¬ 
clusion directly supports the statement that the 
phenomenon of traveling waves in a penumbra is 
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Figure 1. Left panel: V-shaped wakes (white contours) from propagating 3-min waves in the radio¬ 
frequency range (NoRH, 17 GHz). Right panel: EUV narrowband 3-min wakes in the corona (white 
contours). Black contours designate the radio source and coronal loops in the broadband emission. The 
background is the sunspot image in the white light. [Syc/i et al., 2010] 


a chromospheric response to the upward propagat¬ 
ing magnetoacoustic waves generated in the pho¬ 
tosphere. 

For transition zone using a short-baseline ra¬ 
diointerferometer (SBRI) [Zandanov and Uralov, 
1983, 1984] showed the existence of oscillation 
sources, distributed in space with 3-, 5- and 7- 
minute periods. Further, using one-dimensional 
observations at the Siberian Solar Radio Tele¬ 
scope (SSRT, 5.7 GHz) [Zandanov and Sych, 1989; 
Merkulenko et al, 1992], as well as the Nobeyama 
radioheliograph (NoRH, 17 GHz) two-dimensional 
images [Gelfreikh et al, 1999], revealed were the 
radio brightness variations leading not only to a 
brightening of the part of the source coinciding 
with a sunspot umbra, but also to periodic spa¬ 
tial displacements. This indicates at the existence 
of a fine spatial structures inside an active region 
in the form of oscillating small-angular-size sources 
that are not spatially resolved by a radiotelescope. 

Further investigations [Shibasaki, 2001] by us¬ 
ing correlation curves and two-dimensional images 
of radio sources showed that the 3-min oscillations 
steadily exist in sunspots over long periods. Apply¬ 
ing the density and temperature values obtained at 
simultaneous observations with SOHO/SUMMER 
to calculate the gyroresonance emission, found 
was a good agreement with the parameters of the 
detected radio emission. It was concluded that 
brightness variations in sources may be caused by 
density and temperature oscillations when acous¬ 
tic waves pass upward through the third gyroreso¬ 


nance level with the corresponding mean value of 
the magnetic field (~2000 G). In this study, inte¬ 
gral radio images of sources were used, which did 
not allow one to obtain their (sources) fine spatial 
and frequency structure. [Gelfreikh et al., 2004] 
suggested that the periodicity might be explained 
by the existence of resonance structures for acous¬ 
tic and Alfv& waves. The position of the res¬ 
onators, as well as their size, may change within 
rather broad limits and, correspondingly, generate 
different observed oscillation periods. Three basic 
resonator types have been proposed: the resonator 
coinciding with the emission region, the resonator 
outside the emission region, but quite close to the 
radio source, and the resonator of the global solar 
nature, similar to 5-min oscillations. 

[Nindos et al., 2002] used investigation into dif¬ 
ference radio maps obtained with the VLA radiote¬ 
lescope in the two wavelengths, 8.5 and 5 GHz, 
with a high spatial resolution. The emission oscil¬ 
lations from a sunspot are shown to indicate sharp 
variations in brightness with a fast rise and slow 
drop. Despite such variations, they feature spa¬ 
tial, amplitude, and phase stability. The inten¬ 
sity variation spatial distribution is presented like 
fine-structure details of a small angular size at the 
umbral boundary. These structures have a circu¬ 
lar distribution as a patches, and their emission 
arrives from the gyrofrequency second level. 

[Sych and Nakariakov, 2008] studied spatial, 
temporal, and phase dependences of oscillations in 
the sunspot-related radio sources (NoRH, 17 GHz). 
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Figure 2. The oscillating power information extracted from the 304 A data set of the sunspot AR11131 
(08 Dec 2010), shown in the polar coordinates. The dashed lines mark the borders of the umbra and 
penumbra. The unreliable region is cross-hatched, a) The peak power distribution, b) The peak period 
distribution, c) The ID spectral power map as a function of distance and period. The contour in asterisk 
shows the cutoff period, d) The spectral power variance distribution, e) The spectral noise distribution, 
f) The reconstruction of the magnetic field inclination. [Yuan et al., 2014] 


Because the instrument resolution is not sufficient 
for direct spatial resolution of wave, the Pixelized 
Wavelet Filtration method (PWF-analysis) was 
developed. This numerical method is based on the 
signal decomposition for each pixel of the image 
cube into individual spectral components by us¬ 
ing wavelet transform. This allowed one to obtain 
not only oscillation power spatial distribution for 
individual harmonics in the oscillation spectrum 
(narrowband sources), but also to trace their time 
dynamics. The 3-min sunspot oscillations are lo¬ 
calized in the center of microwave sources related 
to sunspot umbra. The 5-min component is local¬ 
ized in the small-size symmetric details at the um- 
bral/penumbral boundary. This location is similar 
to that observed with VLA [Nindos et al, 2002], 
taking into account the difference in the spatial 
resolution of instruments. There are signal phase 


differences depending on the spatial localization of 
details. 

[Sych et al, 2010] found that the 3-min os¬ 
cillation power distribution in the radio sources 
(NoRH, 17 GHz) like a patches inside umbra 
and coincide with the foot points of open mag¬ 
netic loops in the TRACE 171 A ultraviolet range 
(Fig.l). The emergence of these spatial details 
in the morphology of the 3-min wave traces can 
be interpreted as the non-stationarity of the wave 
leakage from the sub-photospheric regions to the 
transition region. The wave traces are found to 
have round and extended shapes. This indicates 
the change of the wave-vector component in the 
plane perpendicular to the LOS. In the quiet stage 
of sunspot activity, the wave traces had a round 
shape, indicating that the waves were propagat¬ 
ing mainly along the LOS or were standing in the 
LOS direction. As the LOS was almost perpen¬ 
dicular to the solar surface in the vicinity of the 
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sunspot, this corresponded to the vertical motions 
of plasma. The periodical increases ~ 10-30 min 
in wave disturbances can result the emergence of 
V-shaped wakes and correspond to the horizontal 
wave propagation in foot points of inclined mag¬ 
netic tubes in the corona. The fine structure of 
the three-minute wave traces at the transition re¬ 
gion and coronal levels is found to almost coin¬ 
cide. Thus, it is likely that these phenomena are 
caused by the same wave motion that propagates 
from the deeper sunspot regions into the corona. 
For 3-min oscillations, there exists a correlation be¬ 
tween the increase in the wave amplitude and the 
increase in the horizontal projected phase speed. 
The change of the speed is likely to be attributed 
to the change in the magnetic-held geometry, per¬ 
haps caused by the waves. In the transition re¬ 
gion the growth of the 3-min amplitude is accom¬ 
panied by the decrease in the LOS angle to the 
wave-propagation direction. In other words, the 
higher amplitude waves are observed to propagate 
more vertically. A possible interpretation of this 
phenomenon is the effect of the centrifugal force 
produced by the field-aligned flows of plasma along 
the curved magnetic-field lines [Nakariakov et al., 
2000; De Moortel, 2006; de Moortel, 2009]. 

For high-frequency oscillations, the sources with 
the oscillation period less than 3 minutes are local¬ 
ized in the umbra and decrease in size as the pe¬ 
riod decreases [Sych and Nakariakov, 2008; Yuan 
et al., 2014; Jess et al., 2012]. As the period grows, 
the size of sources also increase. At the chromo¬ 
sphere level, they reach the maximal value with 
the oscillation period of about 20 minutes, localiz¬ 
ing at the penumbral boundary. For the symmetric 
sunspots, the sources of low-frequency oscillations 
with the period more than ~ 3 minutes look like 
expanding rings, resembling the umbral boundary 
in shape. Small bright patches of the 5-min oscilla¬ 
tion sources in the radio-frequency range (NoRH, 
17 GHz) [Sych and Nakariakov, 2008] and (VLA, 
8.5 GHz) [Nindos et al., 2002] may merge, forming 
circles [Yuan et al., 2014; Kobanov et al., 2013a; 
Zhugzhda et al, 2000]. The localization of the 
sources at the umbral/penumbral boundary indi¬ 
cates at a strong coupling between acoustic waves 
and the magnetic field in the region, where local 
conditions favor the absorption of the global solar 
p-mode [Colly et al., 2003; Schunker and Cally, 
2006]. The observed spatial dependencies of os¬ 
cillation sources can be explained by the spatial 
distribution of the magnetic field and its curvature 
in the framework of modification of the cut-off fre¬ 
quency [Bel and Leroy, 1977]. 

[Reznikova et al, 2012] found that the SDO/AIA 
ultraviolet high-frequency oscillations within the 
1.8-3.3 min range are mostly pronounced closer 


to the umbral center, while the lower-frequency 
component in the form of oscillating rings concen¬ 
trates on peripheral sites of the umbral/penumbral 
boundary. Similar spectral characteristics of oscil¬ 
lations were found for all the temperature chan¬ 
nels. The oscillation region slightly expands with 
height. The oscillating ring diameter broadens by 
about 1 arcsec with each subsequent SDO/AIA 
temperature channel, overstepping the boundary 
at 304A. In the 171A-211A coronal channels, the 
high-power oscillations are localized in coronal fan 
structures. Higher frequencies concentrate toward 
the center at 1.8-2.1 min range, but their power de¬ 
creases. The cut-off frequency is correspondingly 
observed to decrease from 5.7 mHz to 5.0 mHz or 
increase from 2.9 min to 3.3 min in period range. 
This discovery is interpreted as a manifestation for 
the variation in the inclination of the magnetic 
field lines that exit the sunspot at the tempera¬ 
ture minimum level. This confirms the assump¬ 
tions that sunspot strong magnetic fields play the 
role of waveguides for acoustic waves propagating 
upward from the subphotospheric levels into the 
corona. 

[Reznikova and Shibasaki, 2012] calculated the 
magnetic field in potential approximation and 
studied the distribution of the oscillation power in 
the sunspot atmosphere from the SDO/AIA data. 
The authors shown that the propagating waves 
are manifestations of a common oscillatory phe¬ 
nomenon, like the leakage of the solar global p- 
mode from the photosphere [De Pontieu et al., 
2005]. The size of the 3-min oscillation sources 
localized in the umbral region was found to grow 
with height as magnetic field lines expand. There 
is transformation of the oscillation period from 3- 
min in the umbra to 5-min traveling waves in the 
penumbra. This effect may be explained by the 
frequency variation in the sunspot region. To in¬ 
terpret this variation, the authors have to refer to 
the theory of the magnetoacoustic gravity (MAG) 
waves, that is, the acoustic-like slow-mode propa¬ 
gation along magnetic field lines. Due to the pres¬ 
ence of a gravitational field, compressible MAG 
waves can propagate upward only at frequencies 
above the cut-off frequency with effective gravity 
dependence. In regions of low plasma /3, which 
is the case above the sunspot umbra, the effective 
gravity on a particular magnetic field line should 
be decreased by the cosine of the inclination an¬ 
gle 0 of that field line with respect to the direc¬ 
tion of gravity. Therefore, variation of the cut-off 
frequency across the sunspot and observed spatial- 
frequency shape of oscillation sources accordingly, 
depends only on the change of the held inclination 
angle. 

The study [Yuan et al, 2014] has become a fur¬ 
ther investigation into the cut-off frequency ef- 
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feet on the distribution of the oscillation source 
shape in the polar coordinates (Fig. 2). The 
oscillations up to 20 minutes were studied. The 
short-period oscillations were shown to concentrate 
within the umbral boundaries. The long-period 
sources of oscillations are located in the penum- 
bral region as expanding rings. There exists a 
dependence between the oscillation power peak in 
the Fourier spectrum and the source boundary dis¬ 
tance to the sunspot center. The lowest-frequency 
oscillations concentrate at the penumbral bound¬ 
ary. The magnetic field was reconstructed based 
on the MAG wave theory and spatial distribu¬ 
tion of the cut-off frequency. A good agreement 
was found with the parameters of the extrapolated 
magnetic field in potential approximation. The 
obtained values for the inclinations of the recon¬ 
structed fields are mainly higher, than those in the 
case of the extrapolated field. This discrepancy 
assumes a significant effect of other physical pro¬ 
cesses affecting propagation of waves, such as ra¬ 
diation cooling that can the change in cut-off fre¬ 
quency [Centeno et al., 2006, 2009; Felipe et al, 
2010; Roberts, 1983; Ulmschneider et al., 1991; Sri- 
vastava and Dwivedi, 2010a] and MHD mode con¬ 
version [Parchevsky and Kosovichev, 2009; Cally 
and Hansen, 2011; Khomenko and Cally, 2012]. 
The developed method is alternative for the field 
reconstruction, and may be applied to specify the 
boundary and/or initial conditions when using the 
routine methods for the magnetic field extrapola¬ 
tion. 

Proceeding from the cut-off frequency distribu¬ 
tion over a sunspot, one may assume that narrow- 
band sources with a period less than 3 minutes 
will also have a ring-shaped structure like lower- 
frequency components [Reznikova et al, 2012]. In¬ 
side the umbra with the magnetic field vertical 
component, a slow variation in the field spatial in¬ 
clination relative to the solar normal line occurs. 
The cut-off frequency, as a cosine of the field line 
inclination angle, slowly changes. Correspondingly 
to these changes, slow variations in the diameter of 
ring-shaped sources will occur in the 3-min oscil¬ 
lation band with a continuous filling of the umbra. 
For the penumbra with a fast change in the field 
inclination, the rate of the change in the source di¬ 
ameter increases, and the detected rings, without 
superposing on each other, are correspondingly ob¬ 
served. 

[Sych and Nakariakov, 2014] obtained in first 
time the evidence of a fine spatial structure of os¬ 
cillation sources in the 1.5 - 3.5 min oscillation 
band, and compared them to the Fourier spec¬ 
trum harmonics. The authors performed spatial 
spectral decomposition of the umbra in the 3-min 
oscillation band by using the PWF-analysis [Sych 


and Nakariakov, 2008]. The observation duration 
was one hour. Figure 3 presents the obtained re¬ 
sults in the form of narrowband difference maps 
for the detected narrowband harmonics (Fig. 3, 
Left panel ), and their spatial spectral dynamics 
on the time-distance plots (Fig. 3, Right panel ). 
Scanning the narrowband source brightness varia¬ 
tion occurred through the umbral center. Fig. 3, 
Left panel shows that the high-frequency oscilla¬ 
tions with the ~1.5-min period (the beginning of 
the 3-min band) mainly concentrate in the center. 
As the periods grow, the rings form, whose diame¬ 
ter variations create V-shaped diverging branches 
(Fig. 3, Right panel ). A minimal filling (point 
high-frequency source) will be observed at the be¬ 
ginning of the 3-min band. The umbra filling max¬ 
imal value will occur at the end of the band. This 
results in the effect of the frequency rings’ nesting 
in one another, and in forming continuous filling of 
the umbra, which is observed for the 3-min band 
oscillation sources [Reznikova et al., 2012]. The 
oscillation power variation irregularities along the 
branches agree with a fine-structure distribution 
of spectral harmonics in the 3-min oscillation band 
of the integral Fourier spectrum. In the penum¬ 
bra, the inclination of magnetic field lines sharply 
increases, which leads to acceleration of the ring 
diameter increase for the lower-frequency compo¬ 
nents. 

[Yuan et al., 2014] analyzed 5-min oscillations 
in the region of umbral light bridges by using the 
HINODE/SOT instrument. This component is 
shown not to vary dramatically along the bridge, 
which indicates at the oscillation emergence be¬ 
neath. The obtained oscillation characteristics are 
similar to those of the penumbral 5-min traveling 
waves. The 3-min oscillations are shown to be lo¬ 
calized in all the umbral mottles separated by light 
bridges. The opposite sides of bridges oscillate in 
phase, which indicates at a uniform source of 3-min 
oscillations for different mottles. Umbral flashes 
are discovered to be a continuation of the umbral 
oscillations, with a larger amplitude, and with no 
effect on the fundamental oscillation phase. The 
oscillation period is about 3 minutes. 

2.1. Summary 

There is a dependence revealed between the 
spatial-height distribution of the cut-off frequency 
in the sunspot atmosphere and the fine structure of 
the oscillation sources in narrowband images. The 
oscillation sources with periods less than 3 minutes 
are shown to be located in an umbra, decreasing in 
size as the cut-off frequency increases. The 5-min 
oscillation sources are localized around the umbra- 
penumbra boundary. With the period growth, the 
size of the sources increase, reaching the penumbra 
boundaries. This regularity reflects the dynamics 
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Figure 3. Left panel: Running-difference narrowband maps of the umbral oscillations signals obtained 
at 304 A, integrated from 02:30 to 03:30 UT, at the frequencies 11.11 mHz (1.5 min), 8.3 mHz (2.0 
min), 6.7 mHz (2.5 min), and 5.2 mHz (3.2 min). Right panel: Dependence of the narrowband power 
distribution across the umbra (in the horizontal direction via the sunspot center, see the dashed arrow 
in the ~3 min map in Left panel) on the oscillation period. The vertical dashed lines show the umbral 
boundaries. The digits indicate the peaks in the Fourier spectra. [Si/c/i and Nakariakov, 2014] 


of the broadband pulses generated at the subpho- 
tospheric level and propagating in diverging mag¬ 
netic field lines with different inclination to the so¬ 
lar normal and, correspondingly, with a different 
cut-off frequency. In the sunspot atmosphere high 
layers, there form the wakes seen like V-shaped 
brightenings at the coronal loop foot points. The 
detected dependences show that sunspot strong 
magnetic fields play a role of waveguides for the 
acoustic waves propagating into the corona. The 
observational values of the frequencies for different 
layers of the solar atmosphere allow obtaining both 
one-dimensional (in the source plane) and two- 
dimensional (by height) distribution of the field 
line inclination angles and, correspondingly, recon¬ 
structing the magnetic field by using the helioseis¬ 
mology method. 

3. WAVE DYNAMICS IN SUNSPOT 

The chromosphere-corona transition region is a 
very variable zone that is not steady in its height 
position and size. Substance motions are observed 
both up and down. Along with variations in the 
fine spatial localization of sources, [Rendtel et al, 
2003] revealed sign-alternating changes in the form 
of the observed wave frequency drifts in sunspot 
details (plumes) at the SOHO/SUMMER spectro¬ 
graph in the CIV, NeVIII lines, and in the EIT 
ultraviolet radiation (17lA). Most oscillations in 
the transition region and lower corona concentrate 


near 2.7 min, with its amplitude modulation in 
the wave train form, 10-20 minutes (sometimes up 
to 40 minutes) long. Variations in the Doppler 
velocity by the line-of-sight and intensity arise in 
the same period interval (2.3-5.5 min) with a pro¬ 
nounced drift toward short periods. The observed 
frequency drifts of velocity oscillations in the chro¬ 
mosphere as well as in the chromosphere-corona 
transition region may be explained in framework of 
model assuming a parallel stratification of sunspot 
atmosphere with a shrinking resonator size. The 
observability and the line-of-sight component of 
the Doppler velocity depend on the direction of 
the magnetic field. In expanding sunspot magnetic 
field the oscillating volume may change its direc¬ 
tion relative to the line of sight while moving along 
the held lines and thus cause a variation of the ob¬ 
served velocity. Oscillations also can reach regions 
of considerably converging hux tubes and provides 
additional changes of frequency. 

[Centeno et al., 2006, 2009], by using linear 
wave equations with radiation cooling, showed 
good agreement between the observed phase de¬ 
lays and wave amplitude height variations, both 
for sunspots and for pores. A detailed study for 
a larger number of the emission lines emerging at 
different heights is presented in [Felipe et al., 2010]. 
Phase differences of propagating waves and an in¬ 
crease in their power are shown to be intercon¬ 
nected. Therefore, studying the cut-off frequency 
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Figure 4. Wavelet analysis of the 3-min oscillation of the microwave signal on May 3, 2005, from 23:55 
to 03:15 UT. Panel a) shows the amplitudes obtained by calculating the inverse wavelet transform of 
the signal in the 120-240 s period range. Panel b) power distribution of the 3-min oscillation trains in 
the wavelet spectrum. Panel c): time evolution of the 3-min oscillation power. The numbers show the 
numeration of the trains. Panel d) shows the wavelet skeleton representing both global (thick lines) and 
local (thin lines) extrema. In the circles we show the signs of the period drifts (positive, negative or 
without drift). Time is shown in UT, power and amplitudes are in arbitrary units.[Syc/i et al., 2012] 


may also provide information on plasma parame¬ 
ters related to radiative losses. 

[Fludra, 2001], using the SOHO/CDS data, 
showed the existence of significant frequency drifts 
in the 2.8 min oscillations above sunspots in the 
lines of the transition region. The typical train 
duration is about 15 minutes. Different spaced 
sources correspond to different frequencies in the 
spectrum. There are rare instants when the de¬ 
tails with oscillations at several frequencies exist 
simultaneously. 

[Christopoulou et al, 2003], by using wavelet 
spectra, showed the existence of the 3-min oscil¬ 
lation dynamics. There are frequency drifts ob¬ 
served, oscillation modes tend to change, to at¬ 
tenuate or augment over short time intervals. This 
behavior was explained by both the connection be¬ 
tween closely frequency-spaced oscillation modes 
and variations in the physical properties of res¬ 
onators. 

[Sych et al, 2012] showed that the sunspot- 
related 3-min oscillations in radio sources (NoRH, 
17 GHz) have a non-monotonic character. There 
exist low-frequency oscillation trains with an ~8- 
20 min period (~13 minutes being the mean pe¬ 
riod value). The interval between trains is ~13- 
60 minutes. The oscillation trains vary both in 


power and in frequency. The signal variation rel¬ 
ative amplitude is ~3-8%. Pronounced signal fre¬ 
quency drifts are observed throughout the evolu¬ 
tion of individual oscillation trains. The begin¬ 
ning and the end of the drifts agree with the be¬ 
ginning and the end of the amplitude modulation. 
The period changes from 90 to 240 seconds. The 
wavelet analysis showed the existence of both pos¬ 
itive and negative frequency drifts. The frequency 
drift velocity for the 3-min fundamental harmonic 
throughout the low-frequency modulation tends to 
increase with height, reaching 4-5 mHz/hr in the 
photosphere, 5-8 mHz/hr in the chromosphere, and 
11-13 mHz/hr in the corona. Figure 4 presents an 
example of the 3-min oscillation frequency drifts in 
the radio-frequency range. 

There exist significant spatial variations in nar¬ 
rowband 3-min sources throughout the frequency 
drift. New oscillation trains and frequency drifts 
in the transition region are related to the emer¬ 
gence of new oscillation sources in the form of V- 
shaped structures. These changes are noted at all 
the heights of the sunspot atmosphere. When spa¬ 
tially comparing the oscillation sources with ob¬ 
servations in the 171A line (SDO/AIA), obtained 
is their coincidence with the foot points of coro¬ 
nal arches. These structures may be interpreted 
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as the emergence of new waveguides along which 
waves propagate upward. A possible explanation 
for the dynamic properties of the observed waves 
is in the effect of the two factors; dispersion evo¬ 
lution of the wave pulses propagating upward, and 
inhomogeneity of the oscillation power distribution 
over an umbra. The evolution of initial broadband 
pulses results in forming oscillation traces [Rae and 
Roberts, 1982]. Because the sunspot cut-off period 
is about 3 minutes, the observed extended oscilla¬ 
tion sources may be explained as the traces be¬ 
hind the pulses propagating along the magnetic 
field [Fleck and Schmitz, 1991; Kalkofen et al., 
1994; Sutmann and Ulmschneider, 1995; Kuridze 
et al., 2009]. Different magnetic waveguides have 
different physical (temperature, density) and ge¬ 
ometrical (inclination, length, width) parameters. 
Correspondingly to the change in the oscillation 
power and the train emergence, waves propagate 
through different magnetic waveguides forming a 
unique cut-off frequency. The maximum power 
prevalence in a detected tube determines the in¬ 
stantaneous frequency. The power variation and, 
correspondingly, the oscillation frequency are re¬ 
lated to the propagation along different magnetic 
tubes. A sequence of similar variations forms fre¬ 
quency drifts. 

There exists an ambiguous opinion concerning 
the interconnection between the wave processes in 
plumes and coronal fan-structures with the 3-min 
oscillations. [Brynildsen et al., 2004] showed that 
the umbral oscillations are confined to small re¬ 
gions related to foot points of coronal loops. There 
was an assumption about the absence of relation 
between oscillations and plumes. However, [Jess 
et al., 2012] found direct agreement of 3-min oscil¬ 
lations with coronal fans that are anchored in the 
umbra. Their foot points are photospheric umbral 
dots where the oscillation power amplification is 
observed. [Wang et al, 2009] recorded 12- and 25- 
min oscillations of intensity and radial velocity in 
fan like coronal structures over active regions us¬ 
ing the Hinode Extreme-Ultraviolet Imaging Spec¬ 
trometer (EIS) data. These oscillations were iden¬ 
tified as propagating slow magnetoacoustic waves. 

Till present, there is no sufficient information on 
the wavefront dynamics at different levels of the 
solar atmosphere and their stability over long pe¬ 
riods. When observing at a moderate spatial reso¬ 
lution in the 17 GHz range [Gelfreikh et al., 1999], 
the 3-min oscillations were noted to be mainly lo¬ 
calized at the radio source center, and look like 
periodically repeating brightness increases. Fur¬ 
ther, at increasing the temporal and spatial reso¬ 
lutions (TRACE, SDO/AIA), the obtained movie 
and time-distance plots clearly indicate at the dif¬ 
ference between the wave fronts at the bottom 


and at the top of the solar atmosphere. At the 
lower layers of solar atmosphere, there is a quasi- 
spherical propagation of wave fronts with the cen¬ 
ter in sunspot umbrae [Sych et al., 2010]. For the 
coronal level, along with the radial propagation 
along low magnetic loops, there emerge individ¬ 
ual detected directions coinciding with the wave 
propagation along high, open field lines. Their in¬ 
vestigations done by using one-dimensional data 
in the form of scanning the active region with a 
spectrograph slit [Kobanov et al., 2006] allowed 
one to reveal one-dimensional symmetric chevrone- 
type wave fronts on time-distance plots. The front 
structure shows spatial variations with time, which 
indicates at the variation in the wave propagation 
velocity along a sunspot. 

[Sych and Nakariakov, 2014] analyzed the dy¬ 
namics of one- and two-dimensional wave fronts. 
There exist instants on the one-dimensional time- 
distance plots, when a significant departure from 
parallelism between propagating wave fronts is ob¬ 
served. There is direct correlation between ID 
wave form and the oscillating signal level. At the 
oscillation maximum, the waves have a symmet¬ 
ric, spherical shape on either side from the sunspot 
center. At a decrease in the signal level, a decrease 
in symmetry occurs, fractures and spatial displace¬ 
ments with front shape breakup are observed. A 
two-dimensional PWF-analysis [Sych and Nakari¬ 
akov, 2008] of the wavefront trains in the 3-min pe¬ 
riod band showed that, for the maximal power os¬ 
cillation trains, characteristic is the wave propaga¬ 
tion in the form of a spiral with two branches hxed 
in a pulsing point source at the umbral center. A 
counterclockwise rotation of the spiral with a suc¬ 
cessive expansion of the branches up to the shadow 
boundary is observed (Fig. 5). As the power de¬ 
creases, a partial collapse of the branches into indi¬ 
vidual components occurs. Quasi-spherical waves 
form. A similar two-dimensional dynamics is con¬ 
firmed by good correlation with one-dimensional 
time variations on the time-distance plots. One 
should note that the 3-min peak has a certain spec¬ 
tral band. Therefore, taking its hne structure and 
the presence of subharmonics into account, one 
may expect differences both in the shape of nar¬ 
rowband oscillation sources and in their dynamics. 
The analysis of fine frequency structure of the wave 
source shape in the 3-min umbral oscillation band 
showed that the observed quasi-helicity is an as¬ 
semblage of narrowband details with different fre¬ 
quency, spatially separated over the umbral space. 
In the umbral central part, a high-frequency puls¬ 
ing point source with a period of ~1.8 min con¬ 
nected to the lower-frequency arc-wise expanding 
details within the 2.6-3.1 min period range. 

The observed global horizontal oscillations may 
be interpreted within the model based on the fre- 
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Figure 5. Temporal dynamics of the narrowband wave fronts throughout one cycle of 3-min umbral 
oscillations at the transition region level (SDO/AIA 304 A). Black dots indicate front borders. The dash 
circle shows the umbral boundary. [Sych and Nakariakov, 2014] 


quency cut-off of propagating waves. According 
to this model [Zhugzhda, 2008; Botha et ai, 2011; 
Sych et al, 2012], the oscillation frequency de¬ 
creases with decreasing of the magnetic waveguide 
inclination angles relative to the solar normal line. 
In the sunspot center where the field is vertical, 
oscillation sources should have a symmetric shape 
and a maximal value of the cut-off frequency [ Yuan 
et al., 2014] which is observed in reality. The mag¬ 
netic field shows no initial helicity [Reznikova and 
Shibasaki, 2012]. One may assume that the mag¬ 
netic field lines anchored in one place have allo¬ 
cated spatial sectors, where the inclination angles 
will be different in the polar coordinate system. In 
this case, the circular symmetry of wave fronts will 
be violated with the formation of individual seg¬ 
ments in the form of the arc-wise frequency sources 
in the 3-min band. The sequential location of these 
frequency sources from the sunspot center forms 
the observed quasi-helicity. 

3.1. Summary 

The solar sunspot emission dynamics in the 
broad frequency band shows a well-defined low- 
frequency modulation of 3-min oscillations in the 
form of trains. The modulation is non-stationary 
by both the power and the frequency. There are 
frequency drifts over trains. The start and the end 
of the drifts coincide with the start and the end 
of the amplitude modulation. The 3-min oscilla¬ 
tion drift rates were found to grow with height. 


The drift start coincides with the emergence of 
fine-structure oscillating details with the maximal 
power at the foot points of coronal arches. The 
observed spatial-time-frequency peculiarities of the 
oscillations may be treated within two mechanisms 
operating simultaneously. The first is related to 
the dispersion evolution of the subphotospheric 
pulses propagating upward along the sunspot mag¬ 
netic field. The other mechanism is related to a 
spatial split of the propagating pulses into several 
individual magnetic tubes with different physical 
and geometric parameters. For low heights, char¬ 
acteristic is a quasi-spherical propagation of wave- 
fronts near the umbra center. As the height aug¬ 
ments, there occurs the oscillation source broad¬ 
ening, and the emergence of a quasi-spiral wave 
propagation. The cut-off frequency dependence on 
the magnetic field spatial distribution (polar an¬ 
gle) along which there occurs wave propagation, 
and on the emission generation height is assumed 
to be the main factor for the wavefront spiral-like 
form emergence. 

4. SUNSPOT WAVES AND FLARE ENERGY 
RELEASE 

The first ideas of flare initialization in active re¬ 
gions were described in [Norman and Smith, 1978]. 
Therein, there was argumentation that the energy 
release process cannot start throughout the entire 
flare volume simultaneously. At first, the flare is 
localized in a small plasma volume of an active re- 
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gion. Then, the energy release with its dissipation 
starts to capture new volumes of an active region. 
At least, two types of agents may lead to such a dis¬ 
sipative process: electron beams and shockwaves. 
These agents may result in the flare emergence at 
a great distance from the initial position of their 
source, causing simultaneous (sympathetic) flares 
in different active regions [Liu et ai, 2009; Zuc- 
carello et ai, 2009]. 

For the first time, these trigger processes were 
numerically modeled in [Karlicky and Jungwirth, 
1989; Odstrcil and Karlicky, 1997]. [Karlicky and 
Jungwirth, 1989] assumed that electron beams, 
passing through the current layers in the magnetic 
reconnection region, generate Langmuir waves. By 
using a numerical model, the authors studied the 
effect of these electrostatic waves on the plasma 
system. Sufficiently strong Langmuir waves were 
found to be able to generate ion-sound waves 
[Bdrta and Karlicky, 2000]. These waves in¬ 
crease the electrical resistance in the current sys¬ 
tem, which results in the dissipative process start. 
Thereby, electron beams may cause a forced mag¬ 
netic reconnection. 

[Dudik et ai, 2014] showed the effect of the mag¬ 
netic loops’ periodic motions on the emergence of 
the so-called sliding reconnection between loops. 
Analyzed were the microwave emission damping 
oscillations of a flare with a period about 12.5 min¬ 
utes. Similar pulsations were observed in EUV 
in the SDO/AIA 335A high-temperature channel 
[Kim et ai, 2012]. The hot loops observed with 
SOHO/SUMMER were shown to be the oscillation 
source in the form of periodic Doppler shifts. Eor 
the first time, an association between slow magne¬ 
toacoustic waves in the flare microwave emission 
and the hot loop oscillations caused by standing 
waves was shown. 

The latest observation shown that such waves 
take place not only in hot and flaring coronal 
loops, but they are well exist in the comparatively 
cooler loops also. The authors [Srivastava and 
Dwivedi, 2010b] using temporal image data from 
EIS/Hinode, for the first time detect the obser¬ 
vational signature of multiple (first and second) 
harmonics of slow acoustic oscillations in the non¬ 
flaring coronal loop. The investigation suggests the 
simultaneous co-existence of standing oscillations 
and propagating bright blobs. Using 1-D loop hy¬ 
drodynamic model and EIS/Hinode synthetic spec¬ 
tra, the observing possibilities of such oscillations 
were explored [Taroyan and Bradshaw, 2008] with 
prediction that slow acoustic oscillations should be 
observed with the imaging as well as spectral anal¬ 
ysis of real Hinode data. 

[Odstrcil and Karlicky, 1997] studied the flare 
initialization by shockwaves. The authors used a 


2D magnetohydrodynamic model with shockwaves 
propagating through the current sheet. One part 
of the waves passed through the layer, the other re¬ 
flected. Eound was the emergence of wave-caused 
plasma flows near the current sheet. Those flows 
led to the magnetic reconnection emergence. One 
may conclude that a relevant reason leading to the 
reconnection emergence in plasma is not only an 
increase in the electric resistance, but also plasma 
flow emergence. 

The 2005 August 22 long-lived flare observed in 
the radio-frequency range (NoRH, 17 GHz) and in 
X-ray (Ramaty High Energy Solar Spectroscopic 
Imager at 25-50 KeV) was studied in [Reznikova 
and Shibasaki, 2011]. An increase in the signal os¬ 
cillation period from 2.5 to 5 minutes during the 
flare was found. Analyzing the quasi-periodic vari¬ 
ations in the loop length and in the plasma tem¬ 
perature during the flare allowed the authors to 
interpret them as an emergence of the second har¬ 
monic of slow magnetoacoustic waves. This mode 
may be interpreted as a loop response to the initial 
impulse. The mode emergence served a trigger for 
periodic energy release. 

In [Heggland et al., 2009, 2011; Chandra et al, 
2015] the increasing 3-min wave oscillations from 
sunspots before the triggering of the jets was stud¬ 
ied. This effect provide the evidence of wave in¬ 
duced magnetic reconnection. Observed dynamics 
can be interpreted as originating and/or contrac¬ 
tion of new flux tubes, which location near the 
neighbouring magnetic field lines may induce re¬ 
connection and formation of jets. 

An association between periodic initialization of 
reconnection and slow magnetoacoustic waves was 
revealed in [Chen and Priest, 2006]. In this case, 
plasma density periodic disturbances in the recon¬ 
nection region are observed. A similar correlation 
between the 3-min oscillations in the sunspot atmo¬ 
sphere, interpreted as slow magnetoacoustic waves 
and the flare short-period 3-min oscillations has 
been recently established [Sych et al., 2009]. In this 
work the first time, a phenomenological coupling 
between sunspot oscillations and quasi-periodic 
pulsations (QPP) of the flare energy release in the 
active region over the sunspot was shown. Ana¬ 
lyzing the microwave emission recorded with the 
Nobeyama Radioheliograph at 17 GHz showed a 
significant increase in the sunspot 3-min oscilla¬ 
tion train power in the active region prior to the 
flare. The flare curve also showed a modulation 
of the radio flux with a 3-min period. A spa¬ 
tial analysis of the 3-min oscillation distribution 
showed the formation of new V-shaped sources as¬ 
sociated with the foot points of coronal loops in 
the 171A line (SDO/AIA). Visible oscillations in 
the microwave range are propagating waves along 
coronal loops anchored in a sunspot toward the 
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flare source. Sunspot oscillations were suggested to 
appear as a trigger of the flare radio emission. Due 
to the curvature of the magnetic held lines, emerg¬ 
ing centrifugal forces may result in the formation 
of transverse oscillations due to the wave motion 
[Zaitsev and Stepanov, 1989]. These oscillations, 
in turn, appear as a flare trigger due to reconnec¬ 
tion emergence near the local magnetic zero-points 
[Nakariakov et al., 2006]. Using this mechanism 
allowed one to explain both the emergence of in¬ 
dividual short bursts (single reconnections), and 
periodic 3-min modulation (forced cascade recon¬ 
nections due to the wave energy supply) during 
flare evolution. 

A statistical analysis of an increase in the 3- 
min sunspot oscillations before solar flare was per¬ 
formed in [Abramov-Maximov et al., 2011]. The 
parameters of radio emission oscillations (NoRH, 
17 GHz) in active regions, and their effect on the 
flare activity were studied. Illustrated by a large 
number of events, an increase in the 3-min oscilla¬ 
tion power is shown to be characteristic of many 
small flares, which confirm the conclusions in [Sych 
et al., 2009]. A difference in the sunspot oscil¬ 
lations before the flare and after the main peak 
of the energy release was revealed. A correlation 
between the variations in the oscillation spectrum 
and the probability of flare emergence was found. 
The obtained results were interpreted within the 
model for the effect of the waves propagating from 
a sunspot on the flare region. 

A further endeavor to investigate the coupling 
of sunspot waves and flares was the study by [Sych 
et al., 2015], where the authors performed a de¬ 
tailed spatial analysis of the magnetic waveguide 
emergence and evolution of propagating waves. 
The 3-min pulsations in the radio-frequency range 
(NoRH, 17 GHz) with the train period of ~12-20 
minutes were shown to have well-pronounced drifts 
of the oscillation period within ~2-4 minutes. The 
drift value decreases before the flare. This indi¬ 
cates at the formation of the spatially allocated 
magnetic rope (channel), along which slow mag¬ 
netoacoustic waves begin to propagate. The ra¬ 
dio pulse source represents a loopy source, whose 
one foot point is anchored in the umbra. The 
waves propagate from the sunspot toward the flare 
region. For the first time, revealed were obser¬ 
vational data on the transformation of longitudi¬ 
nal low-frequency 3-min waves propagating from a 
sunspot into transverse high-frequency loop oscil¬ 
lations with the ~30-40 sec period (Fig. 6). In the 
coronal region (SDO/AIA, 17lA), the source of 3- 
min EUV oscillations coincides with a radio loop, 
and forms a fine-structure rope of magnetic loops 
(waveguides) connecting the sunspot to the flare 
region. A possible heating mechanism for the flare 


loop is shown to be the coupling (reconnections) 
of the initiated upper-tube transverse oscillations 
with the underlying twisted loops. As the recon¬ 
nection region expands, the magnetic held configu¬ 
ration becomes simpler with the primary loop twist 
disappearance. A flat stage was revealed on the 
flux profile in soft X-ray after the main peak of 
flare emission. Its emergence is assumed to relate 
to maintaining additional energy release due to the 
mechanism for the forced quasi-periodic reconnec¬ 
tions with the underlying loops, initiated by the 
flare loop transverse oscillations. 

4.1. Summary 

There is a relation between the increase in the 
sunspot wave activity and the emergence of flares. 
Revealed was the existence of the sunspot atmo¬ 
sphere wave activity different phases. At the quiet 
phase, there occurs a continuous propagation of 
slow magnetoacoustic waves along the magnetic 
flux tubes (often with an open configuration) from 
the subphotospheric layers into the corona. The os¬ 
cillation sources coincide with the foot points coro¬ 
nal loops anchored in the umbra. For the active 
phase, characteristic is a significant increase in the 
oscillation power within the microwave and an ul¬ 
traviolet radiation with 2-4 min periods. There 
was found the emergence of peak oscillations (by 
power) 20-30 min prior to the beginning of the 
flare energy release. For this phase, characteristic 
is the emergence of V-shaped oscillation sources in 
the umbra. The wave propagation direction coin¬ 
cides with the flare source position. The slow low- 
frequency MHD-waves propagating from sunspots 
were revealed to be able to transform into transver¬ 
sal high-frequency oscillations, and to initiate the 
start of a forced reconnection in magnetic struc¬ 
tures. The waves are shown to possibly be both 
the flare trigger and the emission modulator as the 
intensified wave flux reaches the flare source region 
via the magnetic waveguide. Also revealed was ad¬ 
ditional plasma heating with a production of flat 
patches at the stage of the X-ray flare decay due 
to continuation of the reconnection caused by the 
sunspot wave activity. 

5. MODEL OF SUNSPOT OSCILLATIONS 

Sunspot oscillations were discovered by [Beckers 
and Tallant, 1969]. Unlike 5-min oscillations in 
the quiet atmosphere, the sunspot oscillation pe¬ 
riod appeared to be equal to 3 minutes. How¬ 
ever, soon it became clear that the oscillations, 
both in the quiet photosphere and in sunspots, are 
not only one-frequency oscillations. The spectrum 
of these oscillations comprises many spectral lines. 
Further, the study of the 5-min oscillation spec¬ 
trum led to emergence of a new subdiscipline in 
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Figure 6. Flare loop forced transverse oscillations in the microwave range (NoRH, 17 GHz). Left 
panel: Space-time plot obtained by scanning the flare loop radio brightness cross section. Right panel: 
Variations in the radio loop two-dimensional structure over transverse oscillations (05:57:00-05:57:50 
\JT).[Sych et al., 2015] 


solar physics, helioseismology. This occurred due 
to the creation of the oscillation p-mode theory and 
fast advance in the observational methods. Helio¬ 
seismology and astroseismology, like classical spec¬ 
troscopy of atoms and molecules, allows one to de¬ 
termine many properties of the Sun and stars use 
the oscillation spectrum. 

An important property of sunspot oscillations 
is that they are not global collective oscillations 
of the entire sunspot as a whole. Signatures of 
this were found earlier, but present-day observa¬ 
tions at a high resolution do not leave any doubts 
that the sunspot oscillations are local, i.e. they 
represent an assemblage of many oscillating ele¬ 
ments. Their spectra differ from one another. The 
sunspot oscillation theory should explain, first of 
all, the spectral composition of the oscillations and 
their locality. Currently, sunspot seismology is re¬ 
ferred to as not classical seismology that allows 
one to determine the sunspot atmosphere structure 
by its spectrum, but rather studying the sunspot 
subphotospheric layers through local helioseismol¬ 
ogy with the oscillation p-modes. Sunspot local 
helioseismology is aimed at creating a model for 
the sunspot subphotospheric layers [Moradi et al., 
2010]. In principle, the combination of the lo¬ 
cal seismology methods with the sunspot oscilla¬ 
tion spectrum analysis should help investigate the 
sunspot structure. 

Unfortunately, there has been no sunspot oscilla¬ 
tion model so far that explain the entire oscillation 
spectrum, and not just 3-min oscillations. Within 
the chromospheric resonator model proposed by 
[Zhugzhda and Locans, 1981], there have been en¬ 
deavors to interpret the 3-min oscillation spec¬ 
trum [Settele et al., 2001; Zhugzhda, 2007, 2008]. 
In this theory, the disturbances of the medium 
caused by irregular plasma motions in the subpho¬ 
tospheric layers may generate waves that perme¬ 


ate into the upper atmosphere only at certain fre¬ 
quencies. The sunspot chromosphere is envisioned 
as the Fabry-Perrot filter for propagating slow 
magnetoacoustic waves in strong magnetic fields 
[Zhugzhda et al., 1983]. Recent numerical calcula¬ 
tions [Botha et al., 2011] showed that different pro¬ 
files of plasma temperature and density in the wave 
propagation medium result in different oscillation 
frequency and efficiency of the wave permeation 
into the corona. Therefore, the sunspot atmo¬ 
sphere horizontal inhomogeneity leads to emerg¬ 
ing a fine spectral structure of oscillations, and 
the oscillation power distribution over a sunspot 
[Zhugzhda, 2008; Syeh and Nakariakov, 2014]. Be¬ 
sides, this structuring may be different at different 
heights. 

In order to understand the influence of magnetic 
fields on the propagation properties of waves for¬ 
ward modeling of waves is required. Such calcu¬ 
lations need a model in magnetohydrostatic equi- 
libriumas an initial atmosphere through which to 
propagate oscillations. In [Khomenko and Col¬ 
lados, 2008] developed a method to construct a 
model in equilibrium for a wide range of physical 
parameters. The method combines the advantages 
of self-similar solutions and current-distributed 
models. A set of models is developed by numerical 
integration of magnetohydrostatic equations from 
the subphotospheric to chromospheric layers. 

Further studies [Khomenko et al., 2008] with 
use two-dimensional MHD simulations have shown 
that photospheric 5 minute oscillations can leak 
into the chromosphere inside small-scale vertical 
magnetic flux tubes. The results of numerical ex¬ 
periments are compatible with those inferred from 
simultaneous spectropolarimetric observations of 
the photosphere and chromosphere obtained with 
the Tenerife Infrared Polarimeter (TIP) at 10830A. 
The efficiency of energy exchange by radiation in 
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the solar photosphere can lead to a significant re¬ 
duction of the cutoff frequency with the propaga¬ 
tion of the 5 minute waves vertically into the chro¬ 
mosphere. 

The urgency of sunspot chromosphere seismol¬ 
ogy is related to that building an empirical model 
of the chromosphere through classical methods 
challenges serious difficulties: the models devel¬ 
oped by different authors [Maltby et al, 1986; 
Staude, 1981] substantially differ from one an¬ 
other [Zhugzhda, 2007, 2008]. However, chromo¬ 
sphere seismology within the chromospheric res¬ 
onator model also faced problems, because it failed 
to explain all the complex spectrum of sunspot os¬ 
cillations. Moreover, these studies were based on 
observations of sunspot oscillations [Bogdan and 
Judge, 2006] at a relatively low spectral resolution, 
which, as it appeared later, did not allow one to 
obtain a detailed oscillation spectrum. Only after 
the SDO/AIA launch, it became possible to ob¬ 
tain spectra at a significantly larger spectral res¬ 
olution. The 3-min oscillation spectrum appeared 
to involve much more spectral harmonics, than it 
had been considered earlier [Reznikova et al., 2012; 
Yuan et al., 2014]. 

[Zhugzhda and Sych, 2014] performed a spec¬ 
tral analysis of the data with a high temporal and 
spatial resolution of the ultraviolet oscillations at 
the transition level (SDO/AIA, 304A). The ob¬ 
tained spectrum of 3-min oscillations involves a 
large number of narrowband harmonics. There is 
no detected single 3-min peak. The oscillations 
concentrate in ~1500 km spatial cells with almost 
identical frequencies, without global oscillations of 
the entire sunspot. Oscillation cells emerge in 
those locations, where advancing plasma, moving 
apart the magnetic field, reaches the sunspot pho¬ 
tosphere lower boundary. In these plasma regions, 
there appears an increased temperature gradient 
and a weakened magnetic field. There exists a set 
of hot plasma fluxes with a weak magnetic field 
surrounded by a relatively strong magnetic field. 
It is possible that the oscillation cells are related 
to bright umbral dots. 

The chromospheric resonator is shown not to be 
responsible for the 3-min oscillation spectrum, be¬ 
cause it cannot provide such a number of spectral 
harmonics within the 1.7-3.5 min period band. Nu¬ 
merical calculations [Settele et al., 2001] show that 
the chromospheric resonator leads to the emer¬ 
gence of only several spectral harmonics in this 
frequency band, but not tens of harmonics like it 
occurs in reality. The existence of subphotospheric 
resonators may result in the emergence of the spec¬ 
trum with a large number of spectral lines. This 
is possible in case of a quite depth-extended res¬ 
onator with a sufficiently low main resonance fre¬ 


quency. In this case, numerous spectral harmonics 
may be a result of excitation high harmonics in 
the low-frequency subphotospheric resonator rep¬ 
resenting a system of two coupled resonators. One 
may expect the existence of association between 3- 
minute and low-frequency oscillations that are also 
observed in sunspots [Bakunina et al., 2013]. 

The obtained diversities of oscillation spectra in 
different parts of a sunspot find a natural explana¬ 
tion within the Parker sunspot model, according 
to which the sunspot magnetic field is split into 
individual magnetic tubes in the subphotospheric 
layers. These magnetic tubes cannot be identical. 
The differences between the thin tubes should con¬ 
cern the magnetic field value, temperature gradi¬ 
ents, and density. These parameters determine the 
properties of the subphotospheric resonator, and 
the cut-off frequency and phase velocity of slow 
magnetoacoustic waves depend on them. In other 
words, the pattern of sunspot oscillation distribu¬ 
tion is the reflection of the magnetic field structure 
above a sunspot. 

5.1. Summary 

The investigations into the 3-min sunspot oscil¬ 
lations with a high spatial and temporal resolution 
(SDO/AIA) showed that the model of the chro¬ 
mospheric resonator [Zhugzhda and Locans, 1981] 
is not able to elucidate the observations. Study¬ 
ing the oscillation spectrum does not corroborate 
the existence of the detected peak near the 3-min 
period. There is a set of narrowband spectral har¬ 
monics within the 2-4 minute range. The sunspot 
has no global oscillations as a whole. For each har¬ 
monic, one can find a correspondence to a small 
angle-size umbra site oscillating at the given fre¬ 
quency. To elucidate the locality of 3-min oscil¬ 
lations, there is a suggestion of a slow-wave reso¬ 
nance that originates in the subphotospheric mag¬ 
netic tubes (the Parker model) and encompasses 
all the sunspot atmosphere layers. In this case, the 
observed numerous spectral lines may be a result 
of high harmonics excitation in the low-frequency 
subphotospheric resonator. 
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